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1. Introduction 

We revisit supersymmetric Yang-Mills theories on the torus to study the vacuum state in con- 
nection with the Witten index [jlj]. The torus geometry is crucial to preserve the supersymmetry. 
The index counts the number of quantum states (fermionic states with a negative sign). Due to su- 
persymmetry, states at non-zero energy occur in fermionic and bosonic pairs, and do not contribute 
to the Witten index. The counting can therefore be reduced to the vacuum sector. 

Witten argued that the wave function is constant on the moduli space of flat connections and 
the degeneracy due to the constant abelian gluino components costs no energy [[I]]. Requiring gauge 
invariance (under Weyl reflections) one has r + 1 (= 2 for SU(2)) invariant vacuum states which 
are all bosonic, |0 >, U\Q >, U r \0 >, where U = Yl a =i~KK e ^ ■ Untin 1998 there was a 
puzzle: It didn't agree with the infinite volume gluino condensate calculations for SO(N>7) and 
the exceptional groups (Additional references are contained in[||].) Witten solved that by 
finding an extra disconnected component in the moduli space for SO(7), using a D-brane orientifold 
construction 0]. This means that their are isolated points for SO(7) and SO(8) and for SO(N>9) 
there is an extra component which is like the trivial component for SO(N-7). One therefore finds 
l+r{SO{N)) + l + r(SO{N -1)) =h{SO(N)),ov l+[N/2]+l+[(N-7)/2]=N-2 (h is the dualCoxeter 
number). This is precisely the infinite volume result! The field theory construction works also for 
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Table 1: Contributions to the index for small volumes (A. Keurentjes, JHEP 05, 001(1999)) 



exceptional groups [Q, ^], and is based on classifying triples, Q; = Pexp(/ L A,Jx,), which satisfy 
CliClj = CljClj. Also here one finds that the index in small volumes equals the large volume results. 

But a technical problem occurs when using periodic boundary conditions. This is the break- 
down of the adiabatic approximation in the reduction of the degrees of freedom to those of the 
classical vacuum. For SU(2) the classical vacuum is defined up to a gauge transformation by the 
set of zero-momentum abelian gauge fields, A; = where is constant and f are the Pauli 
matrices. Its gauge invariant parametrization is in terms of the Wilson loops that wind around the 
three compact directions of the torus, which are compact variables. This describes the vacuum val- 
ley as an orbifold, T 3 /Z2 for SU(2). Here = exp(— 2nin -xXj/L) (which maps c? to c? + 47T«,) 
and g = 1X2 (which maps c? to — c?) give Gribov copies. The orbifold singularities arise where the 
flat connection is invariant under (part of) the Weyl group (the remnant gauge transformations that 
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leave the set of zero-momentum abelian gauge fields invariant). For SU(2) their are eight orbifold 
singularities (related to A = by anti-periodic gauge transformations), c? = — cf -\-4nrij, giving 
c? = 2ntij (with tij = or 1). It was studied m^fy and we will review some of it here. 



2. The Hamiltonian 

We choose the dependence on the bare coupling constant such that 

^ = -^ 1 (F^ v ) 2 + ^ 1 l a UDnr. (2.1) 

The reduction to the zero-momentum degrees of freedom, as in the bosonic case, will replace the 
bare coupling constant by a running and asymptotically free coupling constant g(L). The zero- 
momentum gauge fields are parametrized as A,- = icfT a /(2L). The vacuum valley is parametrized 
by the abelian degrees of freedom. These are defined by r,-, with r,r y - = £ a c?c", for each i and j. 
As said, we may also parametrize the vacuum valley by r ( - = C, = c]. 



m i > 
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Figure 1: A two dimensional slice of the vacuum valley along the (C\,C2) plane. The grey square is the 
fundamental domain A. The dots are gauge copies of the origin (which turn out to lie on the Gribov horizon 
£1, indicated by the fat square). 

The cell C G [— n,7l] 3 can be used as a fundamental domain A. Any point on the vacuum 
valley can be reached by applying suitable gauge transformations. This reduces the eight orbifold 
singularities to one, but at a cost: Opposite sides on it's boundary are identified by homotopically 
non-trivial gauge transformations (indicated by g in Fig. [[]). The representations of their homotopy 
define the electric flux quantum numbers as introduced by 't Hooft[Q]. We will here only consider 
the sector with zero electric flux, i.e. the trivial representation, where wave functions at opposite 
sides are equal. Since the orbifold singularity implies massless excitations for the transverse fluc- 
tuations, we seemingly have no choice: We have to keep the non-abelian zero-momentum modes, 
and thus the adiabatic approximation breaks down. 

The energy gap in the fluctuations transverse to the vacuum valley is easily read off. Close to 
the origin it is given by 2\C\/L. Integrating out transverse degrees of freedom is only reliable if the 
energy of the low-lying states is smaller than this gap. This energy behaves as g 2 ' 3 (L)/L. Consider 
now a sphere of radius g 1 / 3 (L) around each orbifold, beyond which the adiabatic approximation is 
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accurate. There, the wave function can be reduced to the vacuum valley and, assuming indeed it has 
zero energy, it will be constant. As long as g(L) is small, we may assume the wave function in the 
neighborhood of the orbifold singularity to respect spherical symmetry. In the bosonic case, once 
the wave function spreads out over the vacuum valley, any spherical symmetry is quickly lost. In 
the supersymmetric case, however, the reduced wave function on the vacuum valley will become 
constant up to exponential corrections at distances much greater than g 2 / 3 (L) from the orbifold 
singularities, which are separated over a distance 271. Instead of insisting that the groundstate wave 
function is normalizable, we should rather insist on its projection to the vacuum valley to become 
constant. Spherical symmetry near the orbifold singularities will dramatically simplify the analysis. 

We will now set up the zero-momentum supersymmetric Hamiltonian, and its reduction to the 
gauge invariant, spherically symmetric sector. This is not new [Q], but we will be able to push it to 
the point where we can explicitly construct a complete basis of states that respect these symmetries. 
We start from the supercharge operators, 

ea=CT «^("W"^) ! fi * = ^S*(-^ +ffl /)- (2 - 2) 

with V" = c"/(g(L)L) and a 7 = T ; (and a the unit) as 2 x 2 matrices. Restricting to the zero- 
momentum modes, both the Weyl spinors X% and Xj* are constant. Lowering indices is done with 
e a p = e d p = —iT2, 8 a b and r]^ = diag (1,1,1,-1) (or 6y) for respectively the spinor, group and 
space-time (or space) indices, and raising of indices is done with the inverse of these matrices. 
Repeated indices are assumed to be summed over, but to keep notations transparent we will not 
always balance the positions of the gauge and space indices. Finally, for zero-momentum gauge 
fields one has Bf = -! l gSi jk E abc Vj'V£. 

In the Hamiltonian formulation the anti-commutation relations {X aa , X h $} = 0, {X a& ,X b ^} = 
0, {X aa ,l b P} = d$ a S ah , with db the unit 2x2 matrix (one has (ff^)* a = e^e a ^(o^)^), give 
{Qa, Qa} = 2{a )aaJ^-2(a i )aaV i a %, where 

% = igEabc (vj^b - ~^ b W\ (2.3) 
is the generator of infinitesimal gauge transformations, and is the Hamiltonian density 

* = -^d\^ + ' m ~ i§ea h A a oa h v;\ (2.4) 

Splitting the Hamiltonian, J d 3 xJ$? = g 2 ^(L)H /L, in its bosonic and fermionic pieces, H = Hg + 
H f , we find with cf = g 2/3 {L)cf and W a = - l 1 s ijk £ ahd c b j c d k that H B = -I (d/dcf) 2 + ± {B^f and 
H f = -ie abd l a a i X b cf. 

The orbifold singularities, other than at c = 0, lie at a distance b = 2ng~ 2 / 3 (L) in these new 
variables c (measured along the vacuum valley where B vanishes). We want to solve for the ground- 
state wave function such that for |c| 3> 1 it becomes a constant, after projecting to the vacuum 
valley. As this boundary condition is compatible with spherical symmetry, i.e. it goes to the same 
constant for all directions on the vacuum valley, we will restrict ourselves to wave functions 
that are spherically symmetric and gauge invariant. We stress again that this is an accidental spher- 
ical symmetry, that holds in sufficiently small volumes. 
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3. Vacuum valley and boundary condition 

States with odd fermion number F do not respect the symmetry and particle-hole duality re- 
lates F = to F = 6 and F = 2 to F = 4 Without fermions (F = 0) this will reduce to the 
bosonic Hamiltonian which is known to have a non-zero vacuum energy^]. We therefore write 
down the most general F = 2 states, that is the two-spinor states which are symmetric or antisym- 
metric in the gauge index (and thus respectively antisymmetric and symmetric in the spinor index), 
\V) = Vj a ^ 3 a = -2i^e ahc l^)fX h ^\0) and \3>) = y ab J? ab = -^ ab lfX b e^\0) (where 
a' JV = \{6^o v — a v a^), such that a ;0 = \%j as a 2 x 2 matrix with the first index up and the 
second down). Here if ? and <¥" tb = 5? are arbitrary and covariance allows us to write this as 
Yf =h i c a j /r-h 2 B a j /r 2 + h 3 c b j c b k c a k /r 3 and J^ ab = h 4 8 ab + h 5 c a jC b /r 2 + h 6 c a j cjcfc b k /r 4 , where h t 
depend on f 2 = {c") 2 , u = r (B") and v = f~ 3 detc. This gives a six dimensional matrix equation 
for the Hamiltonian which can be split into a "radial" and "angular" part (Hf = Hf/r), 

H^ = (p\n'\H\p,n)=E e p »8^8PP' + (p\£ n ,\f*\p,£„)tf (3.1) 

Our basis explicitly diagonalizes the angular part ("spherical harmonics") in terms of invariant 
polynomials. In Fig. § the angular band structure is given by plotting (n'\u\n) ^ or (n'\Hf\n) / 




Figure 2: Band structure of the reduced Hamiltonian, for the first 100 states. 

as black squares. The different bands can be traced to come from the selection rules \8£\ = 0,2,4 
for the matrix elements of u, and \8£\ = 1 for the matrix elements of Hf. 

The radial wave function is <j>p(r) = C e p f~ 3 je(Jc p f), with E e p = j(k e p ) 2 and i = 2L + 3 (Je(z) is 
the spherical Bessel function and C p is the normalization constant). It diagonalizes the kinetic term, 
-lr- 8 d f r 8 d f + r~ 2 L(2L + 7) = -^r~ 3 [r- 2 d f f 2 d f + ±r~ 2 ^+ 1)] r 3 . To complete the construction 
of the basis, we need to address the question of boundary conditions (fixing the momenta k? p ), so 
that the groundstate wave function becomes a constant, after projecting to the vacuum valley. The 
vacuum valley is characterized by those configurations for which u = (this implies v = 0), and r 
measures the (three-dimensional) distance to the origin along this vacuum valley. 
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The wave function can be decomposed as *F = £™ =0 r -3 'f n {r)Xn(u,v;r), with #„(w,v;r) nor- 
malized eigenf unctions, so that = Y,7=oJ ^ n ^d?j%{r)f n (f). This means that we have to 
impose d?fo(r) = at the boundary of the fundamental domain, r = b = 7ig~ 2 / 3 (L). This is equiv- 
alent to (xo\dp(f 3x ¥)) = 0, where the inner product is at fixed f. That this boundary condition 
receives corrections is however an artifact of the truncation to the zero-momentum modes. If 
we take into account that in the full theory Xo also involves the non-zero momentum modes, the 
(gauge) symmetry guarantees^] that at the boundary of the fundamental domain d?fo(b) = 0, 
and this source of the breaking of supersymmetry is absent. Higher order terms in computing the 
effective Hamiltonian are required to deal with this [^p. We therefore impose df(f 3v F(f)) = at 
f = b = ng~ 2 ' 3 (L), knowing that it is only valid at small g(L). 

A direct measure for the failure of the Born-Oppenheimer approximation is given by f 2 (r) — 
/q (r) = E" = i/^(?), which as expected, deviates from zero when the gap is small, but also when 
we approach the boundary at f = b (due to the zero-momentum approximation). In Fig. || we show 
the results for b = 4.4, 4.7 and 5.0, to illustrate that this second deviation decreases with increasing 
b (or decreasing coupling). This shows that the mismatch between the boundary condition and 
the truncation of the effective Hamiltonian does not affect the wave function in the neighborhood 
of the orbifold singularities, where the failure of the adiabatic approximation is non-perturbative. 




1 2 3 4 5 2.5 3 3.5 4 4.5 5 



Figure 3: On the right is shown f 2 {f) — /q (f) extracted from the groundstate wave function ^o, satisfying 
the boundary condition df (r iy ¥o)(b) — 0, for b — 4.4, 4,7 and 5.0. We normalized with respect to b — 5, such 
that at f — 2 all f 2 agree. On the left is shown the groundstate energy itself, using up to 20 radial modes for 
each of the 420 harmonics together with the lower bound from Temple's inequality (indicated by the dots). 

Finally Fig. || also shows the groundstate energy Eq (the curve gives the upper and the dots the 
lower bound). To keep the size of the matrix manageable, the components of the eigenvectors 
are removed when they are in absolute value below a threshold, without significantly affecting the 
accuracy. This process of pruning is performed iteratively, increasing the number of radial modes 
per angular state. Together with Temple's inequality [ |l0| ] this is extremely efficient to optimize the 
accuracy and achieve numerical control. 

We have only considered SU(2) to illustrate how to go beyond the adiabatic approxima- 
tion. Our analysis shows that, although the orbifold singularities are cause for concern, in the 
end they do not upset the result for the Witten index. For other groups this is much harder. 
Of course, a numerical analysis can never be entirely conclusive in deciding a theoretical issue 
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that involves the counting of exact zero-energy states. The computer code is still available at 



www.lorentz.leidenuniv.nl/vanbaal/susyYM (partly they could be used for other studies Qllfl). 
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